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Abstract 

We address the problem of the response to a seismic wave of an ur- 
ban site consisting of a large and infinite number (iV&) of identical, eq- 
uispaced blocks overlying a soft layer underlain by a hard substratum. 
The results of the theoretical analysis, appealing to a space- frequency 
mode-matching (MM) technique, are compared to those obtained by a 
space-time finite element (FE) method. The two methods are shown to 
give rise to much the same prediction of seismic response for Nb = 10. 
The MM technique is also applied to the case Nb = oo, notably to 
reveal the structure and natural frequencies of the vibration modes of 
the urban site. The mechanism of the interaction between blocks and 
the ground, as well as that of the collective effects of the blocks, are 
studied. It is shown that the presence of a large number of blocks 
modifies the seismic disturbance in a manner which evokes, and may 
partially account for, what was observed during many earthquakes in 
Mexico City. Disturbances at a much smaller level, induced by a small 
number of blocks are studied in the companion paper. 
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1 Introduction 



The Michoacan earthquake that struck Mexico City in 1985 presented 
some particular characteristics which have since been encountered at the 
same site and at various other urban sites [SH [H3 \T7\ H21 > but usually at a 
smaller scale. Other than the fact that the response in downtown Mexico 
varied considerably in a spatial sense ^1 EH] , was quite intense and of very 
long duration at certain locations (as much as «3min |30l 1331 I15j). and 
often took the form of a quasi-monochromatic signal with beatings |28| I33j . 
a remarkable feature of this earthquake (studied in El E3 ) was that 
such strong motion could be caused by a seismic source so far from the 
city (the epicenter was located in the subduction zone off the Pacific coast, 
approximately 350km from Mexico City). It is now recog nized E] that 
the characteristics of the abnormal response recorded in downtown Mexico 
were partially present in the waves entering into the city (notably 60km 
from the city as recorded by the authors of ^S]) after having accomplished 
their voyage from the source, this being thought to be due to the excitation 
of Love and generalized-Rayleigh modes by the irregularities of the crust 

0G3EES]). 

In the present investigation (as well as in the companion paper), we 
focus on the the presence of the built features of the urban site as a com- 
plementary explanation of the abnormal response. In the companion paper, 
we treat the case of one or two (different or identical) built features in the 
form of cylindrical blocks. Herein, we study the case of many (10, 20, 40,...., 
oo) identical blocks. Such a configuration is a more realistic representation 
of a real city due to the large number of blocks it incorporates, but the 
assumptions that the blocks are: i) cylindrical, ii) identical and iii) peri- 
odically arranged, are rather far removed from reality, except in restricted 
portions of modern cities and megacities. Nevertheless, these assumptions 
are not more unrealistic than random dispositions of blocks with random 
sizes I34| I17j and compositions, and have the advantage of enabling a 
theoretical analysis which can shed some light on the physical origins of the 
above-mentioned exotic phenomena. 

The periodic model of cities built on sites with a soft layer overlying a 
hard substratum solicited by seismic waves originated in the work of Wir- 
gin and Bard |38j . Unfortunately, the theoretical apparatus underlying the 
numerical results was not given in this paper (only references were made 
to studies that treat this issue, but as most of these references were rela- 
tive to problems of electromagnetic waves, they have escaped the attention 
of the civil engineering and geophysical communities) and the distance be- 
tween blocks in the computational results was taken to be 2000m, which, 
in our present opinion, is unrealistically large (unless the blocks represent 
skyscrapers, of which there are usually few, and largely-distant one from the 
other). 
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The object of the present investigation is thus twofold: i) give the missing 
theoretical foundations of the results of _3Sj , and ii) treat cities that are more 
realistic than those in |38j and in the companion paper. More specifically, 
we shall address the following questions: 

(i) how should one account for the principal features of the seismic re- 
sponse in the cases of a relatively-large number of blocks (the case of a 
relatively small number of blocks being treated in the companion paper)? 

(ii) what are the vibrational modes of the global structures (i.e. the 
superstructure plus the geophysical structure) and what are the mechanisms 
of their excitation and interaction? 

(iii) what are the repercussions of resonant phenomena on the seismic 
response? 

(iv) what are the differences in seismic response between configurations 
with a small and a large number of blocks? 
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Figure 1: Satellite picture of a portion of Mexico City, Mexico. 




Figure 2: Satellite picture of a portion of Nice, France. 



2 Candidate sites 

Many earthquake-prone cities and megacities (New Delhi, Tokyo, Mexico 
City, Istanbul, San Francisco, Basel, etc.) are built on soft soil underlain 
by a hard substratum and contain districts with periodic, or nearly-periodic 
arrangements of blocks or buildings. An attempt will be made in the present 
investigation to analyze the seismic response of districts of this type in the 
two cities, Mexico City and Nice (satellite pictures of which are depicted in 
figs, n and El respectively) . 
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Figure 3: View of the 2D city (only two of the blocks are represented). 

3 Description of the configurations 

We focus on a portion of a city characterized by a periodic or quasi- 
periodic set of identical blocks, assumed to have 2D geometry, with X3 the 
ignorable coordinate of a Ox\Xix?, cartesian coordinate system (see fig. IB}. 
The blocks are in welded contact with the ground underneath which is lo- 
cated a horizontal soft layer underlain by a hard half space (see fig. OJ) • Each 
block is composed of one or more buildings, all the blocks have the same 
shape (rectangular in the cross-section plane), are of the same size (charac- 
terized by two constants, their height b and width w) and composition, and 
their separation (center-to-center spacing) is a constant d. For the purpose 
of the analysis, each block is homogenized, so that the final geometry of the 
city is as represented in fig. 0J Fig. El represents a cross-section (sagittal 
plane) view of the city, is the stress- free surface composed of a ground 
portion T g , assumed to be flat and horizontal, and a portion T ag , consti- 
tuting the reunion of the above-ground-level boundaries Tag ; j £ 7L of the 
blocks. The ground Tq is flat and horizontal, and is the reunion of T g and 
the base segments Ti ; j £ Z joining the blocks to the underground. 

The medium in contact with, and above, is air, assumed to be the 
vacumn (i.e., Tf is stress- free) . The medium in contact with, and below Tq 
is the mechanically-soft layer occupying the domain fij, which is laterally- 
infinite and of thickness h, and whose lower boundary is Th, assumed to 
be flat and horizontal. The soft material in the layer is in welded contact 
across with the mechanically-hard material in the semi-infinite domain 
(substratum) Qq. 

The domain of the j-th block is denoted by fig an d the reunion of all the 
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Figure 4: View of the 2D city with homogenized blocks (only two of the 
blocks are represented). 




Figure 5: Sagittal plane view of the 2D city with a set of identically, equally- 
spaced, homogenized blocks. 
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^2 ; j G B is denoted by f^. The material in each block is in welded contact 
with the material in the soft layer across the base segments T 3 bs ; j'gZ. 

The origin O of the cartesian coordinate system is on the ground, 22 
increases with depth and 23 is perpendicular to the (sagittal) plane in figs. 
13151 With ij the unit vector along the positive Xj axis, we note that the unit 
vectors normal to Tq and I\ are — 12. 

The media filling Q , Q 1 and Sl{ are M°, M 1 and M 2 i = M 2 ; V j G B 
respectively and the latter are assumed to be initially stress-free, linear, 
isotropic and homogeneous (thus, each block, which is generally inhomoge- 
neous, is assumed to be homogenized in our analysis). We assume that M° 
is non-dissipative whereas M 1 and M 2 are dissipative. 

The seismic disturbance is delivered to the site in the form of a shear- 
horizontal (SH) plane pulse wave, initially propagating in Qq. The SH nature 
of the incident wave (which fact is indicated by the superscript % in the 
following) means that the motion associated with it is strictly transverse 
(i.e., in the 23 direction and independent of the 23 coordinate). Both the 
SH polarization and the invariance of the incident wave with respect to 23 
are communicated to the fields that are generated at the site in response to 
the incident wave. Thus, our analysis will deal only with the propagation 
of 2D SH waves (i.e., waves that depend exclusively on the two cartesian 
coordinates x%, 22 and that are associated with motion in the 23 direction 
only). 

We shall be concerned with a description of the elastodynamic wavefield 
on the free surface (i.e., on Tf) resulting from plane seismic wave solicitation 
of the site. 

4 Governing equations 

4.1 Space-time framework wave equations 

In a generally-inhomogeneous, isotropic elastic or viscoelastic medium 
M occupying M 3 , the space-time framework wave equation for SH waves is: 

V • ( M (x, w)Vu(x, t)) - p(x)^ 2 n(x, t) = -p(x)/(x, t) , (1) 

wherein u is the displacement component in the i3 direction, / the compo- 
nent of applied force density in the i3 direction, \x the Lame descriptor of 
rigidity, p the mass density, t the time variable, u the angular frequency, df 
the n— th partial derivative with respect to t, and x = (21,22). Since our 
configuration involves three homogeneous media and the solicitation takes 
the form of a plane wave initially propagating in Qq (thus, the source density 
is nil for this type of wave), we have 

(c"») 2 V • Vu m (x, t) - dfu m (x, t) = ; x G O m ; m = 0, 1, 2j , j G B , 

(2) 
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wherein m superscripts designate the medium (0 for M°, etc.), and c m is 
the generally- complex velocity of shear body waves in M m , related to the 
density and rigidity by 

(c m (-)) 2 = ^ , (3) 

it being understood that p m , fj> m (uj) ; m = 0, 1, 2j, j £ B are constants with 
respect to x. In addition, the densities are positive real, and we assume 
that the substratum is a dissipation-free solid so that the rigidity therein is 
a positive real constant with respect to to, i.e., /^ (w) = p° > 0. 

4.2 Space-frequency framework wave equations 

The space-frequency framework versions of the wave equations (j2J are 
obtained by expanding the displacements in Fourier integrals: 

/oo 
u m (x,uj)e~ lujt du) ,Vt € R , (4) 
-oo 

(wherein i = \/— I, o> the angular frequency and i the time) so as to give 
rise to the Helmholtz equations 

V-Vn m (x, uj)+(k m (uj)) 2 u m (x, w) = ; Vx£fi m ; m = 0, 1 , 2j, j e B , 

(5) 

wherein 



k m (uj) := — ^— =UJ -P— . (6) 

v ; c m (a;) y /i m (u;) v ; 

is the generally-complex wavenumber in M m . 

We shall deal with constant quality factors Q m ; m = 1,2 over the 
frequency range of solicitation, so that the frequency-dependent rigidities 
take the form given in 26 , 



//^H/C/f — ) ;m = l,2, (7) 

wherein: oj re f is a reference angular frequency, chosen herein to be equal to 
9 x lCT 2 Hz. Hence 

/-iuA ^ arctan (o^) 
c m H=C/ ;m = l,2, (8) 



with c™^ := y-pf- Even though Q m , m = 1,2 are non-dispersive (i.e., 
do not depend on w) under the present assumption, the phase velocities 
c m ; m = 1,2 ore dispersive. 
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Due to the assumptions made in sect. 01 



fcV):=£= W ^, (9) 
(i.e., Q° = oo so that k° is real). 

4.3 Space-frequency and space-time framework expressions 
of the incident plane wave 

As mentioned above, we shall be concerned with plane wave excitation 
of the city. 

Actually, a plane wave satisfies the homogeneous wave equation in the 
space-time framework and a homogeneous Helmholtz equation in the space- 
frequency framework. 

The field is chosen to take the form of a pseudo Ricker-type pulse in 
the space-time framework, whose shape is directly connected with the site 
we will consider in the numerical application (either a Nice-like site, or a 
Mexico city -like site) . 

4.3.1 Space- frequency- framework representation of the plane, im- 
pulsive, incident wave 

The plane wave nature of the incident wave is embodied in the choice 

«*(x, u) = S(u) exp (xis* - x 2 d)] , (10) 

wherein s 1 = sin# l , c l = cos8 l and 9 l is the angle of incidence in the x\ — x 2 
plane with respect to the x 2 axis. 

The fact that the incident wave is a pseudo Ricker pulse means that the 
amplitude spectrum S(uj) is given by 



to which corresponds the temporal variation (Fourier inverse of S(u)): 

S(t) = -[2a 2 (t s -t) 2 -l]exp[-a 2 (t s -t) 2 ] , (12) 

for Nice-like site solicitation and by 

ru , 2a 2 uj 2 ( lo 2 \ 

S ^ = -^4^^{ **>-Wj ' (13) 

to which corresponds the temporal variation (Fourier inverse of S(u>)): 

S(t) = 2a 2 [2a 2 (t s - t) 2 - l] exp [-a 2 (t s - t) 2 ] , (14) 
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for the Mexico-like site solicitation. 

In both cases, a = n/t p , t p is the characteristic period of the pulse, and 
t s the time at which the pulse attains its maximal value. In particular we 
will chose t p = t s = where vq is the central frequency (in Hz) of the 
spectrum of the pulse. 

4.4 Boundary and radiation conditions in the space-time 
framework 

Since our finite element method I17j for solving the wave equation 
in a heterogeneous medium M (in our case, involving three homogeneous 
components, M°, M 1 and M 2 ) relies on the assumption that M is a contin- 
uum, it does not appeal to any boundary conditions except on Tf where the 
vanishing traction condition is invoked. The latter is modeled with the help 
of the fictitious domain method which allows us to model diffraction 
of waves by a boundary of complicated geometry, not necessarily matching 
the volumic mesh. Furthermore, since the essentially unbounded nature of 
the geometry of the city cannot be implemented numerically, we take the 
geometry to be finite and surround it (except on the Tf portion) by a PML 
perfectly-matched layer ^2] which enables closure of the computational do- 
main without generating unphysical reflected waves (from the PML layer). 
In a sense, this replaces the radiation condition of the unbounded domain. 

4.5 Boundary and radiation conditions in the space-frequency 
domain 

The translation of the stress-free (i.e., vanishing traction) nature of Tf = 
F g U r ag, with T ag := UjeB r a 9 , is: 

AiVJ^yfow) = ; xer 9 , (15) 

» 2 (u)d n u 2j { X , w) = ; x € r{ g , j € B , (16) 

wherein n denotes the generic unit vector normal to a boundary and d n 
designates the operator d n = n • V. 

Since M 1 and M 2 are in welded contact across := UjgB^L» ^ ne 
displacement and traction are continuous across T^: 

u\*.,u) - u 2 \^uj) = o ; xerj slJ ei, (17) 

nH^dnuH^oj) - fx 2 (uj)d n u 2 3{x,uj) = 0;xGrj slJ £l. (18) 

Since M 1 and M° are in welded contact across I\, the displacement and 
traction are continuous across this interface: 

u 1 (x,w)-u (x,w) ; xer ft , (19) 
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^{^dn^i^u) -/(u;)«9 n u (x,u/) ; x G r fc . (20) 

The uniqueness of the solution to the forward-scattering problem is assured 
by the radiation condition in the substratum: 

it°(x,u;) — u l (x,o>) ~ outgoing waves ; ||x|| — > oo, X2 > h . (21) 

4.6 Recovery of the space-frequency displacements from the 
space-time displacements 

The spectra of the displacements are obtained from the time records of 
the displacements by Fourier inversion, i.e., 

1 r°° 

u m (x, u) = — / u m (x, ije^di ;j = 0, 1, 2j, jGB. (22) 

5 Field representations in the space-frequency frame- 
work for Nb = oo 

Since the formulation for the case of Nf, < oo identical (or non-identical) 
blocks was given in the companion paper, it will not be repeated here. 

The new feature in the present investigation is the possibility of the exis- 
tence of an infinite number (i.e., iVj, = oo) of identical (in shape, dimensions 
and composition) blocks, separated (horizontally) one from the other by the 
constant distance d (called the period). 

At present, the blocks are identified by indices in the set B = {..., — 1, 
0, 1, ...} = Z, with the understanding that the center of the segment T® s is 
at the origin. 

Owing to the plane wave nature of the incident wave, the periodic nature 
of Tq, and the fact that the blocks are assumed to be identical in height, 
width, and composition, one can show that the field is quasi-periodic (this 
constituting the so-called Floquet relation, i.e., 

u j ( Xl + nd, x 2 ) = u j ( Xl ,x 2 )e ik > d ; Vx G % ; Vn G Z ; j = 0, 1 . (23) 
wherein k\ := k°s\ 

5.1 Field in Q 

By referring to the companion paper, and after use of the Green's second 
identity, the field in f^o can be shown to take the form: 

/oo 
B°{ki,u) exp {i \k\X\ + k 2 {x2 — h)] } dk\ ; 
-oo 

xe^o, (24) 
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wherein: 



B (fc 1 ,a;) = - i - i;ff / {d yi u ( yi ,h,u)+ik%u ( yi ,h,u;)}x 

exp (-ik 1 y 1 )dy 1 , (25) 



By a suitable change of variables and use of the Floquet relation B°(ki) can 
be cast in the form: 

i f d/2 

B°(h) = — J / {^ /1 u ( J /i ) fc,a>) + ifc§« (j/i J /i, £ «;)} x 

2 " — d/2 

oo 

exp [— ifciyi] rfyi exp \ixnd (k\ — /ci)] . (26) 

n=— oo 

With the help of the Poisson summation formula it follows that: 

oo oo 

exp[ind{k\-k 1 )]=— Hhn-h) , (27) 



n=— oo 



wherein <5( ) is the Dirac delta distribution and 

fcl„ = k{ + — . (28) 

On the other hand, F(x)5(x — y) = F(y)5(x — y) in the sense of distributions, 
so that 

oo 

= E KHhn-h) , (29) 

n=— oo 

where 

i /" d/2 

^ = ^n / {$aU°(i/i>M) +i^ n u°(yi,/i,w)} x 

iK 2n a J-d/2 

exp (— ifci^yi) dyi , (30) 

and 

*4 = vW-*?„ ; > > > , w > ; i = 0,1 . (31) 

The introduction of (|29|) into ()24() and the use of the sifting property of the 
Dirac delta distribution result in: 



u° 



(x,uj) = u 1 (x,uj) + ^ -B° ex P {i [hnXi + h% n (x 2 - h)] } ; x G Sl • 

n=— oo 

(32) 



14 



However, we can write 

oo 

^(x, oo)= A i ex P I 1 hnxi - k° 2n x 2 ] } ; x£R 2 , (33) 

n=— oo 

wherein 

A\ = S(u;)5 n0 , (34) 

5 nm is the Kronecker symbol, fcio = k\, and A;^ = yj (kj) 2 — (k\) 2 := k^, so 
that 



m°(x,cj)= ^ Al l ex.p{i[k ln x 1 -k% n x 2 ]} + 

n=— oo 

oo 

5° exp {i [fe ln s! + k° 2n {x 2 - h)] } ; x G fi , (35) 

n=— oo 

with the understanding that a := {A^ ; n G Z} is a known vector and 
b° := {2?° ; n £ Z} an unknown vector. 

5.2 Field in fii 

By proceeding in the same manner as previously, we find: 



u 



1 (x,w)= A^exp{i[/ci n xi-^ n x 2 ]} + 

n=— oo 

oo 

^ B\ exp {i [fei„xi + k\ n x 2 ] } ; xefli, (36) 

n=— oo 

with the understanding that both a 1 := {A* ; n G Z} and b 1 := ; n G 
Z} are unknown vectors. 

5.3 Field in Q 2 

The Floquet relation actually extends to Vt 2 wherein it takes the form 

u 2j (x 1 ,x 2 ,io) = u 20 ( Xl ,x 2 ,u;)e ik i jd ; j G Z . (37) 

By referring to the companion paper, the displacement field in the zeroth-th 
block takes the form: 

oo 

m 20 (x, u) := u 2 (x, u) = B l cos \k\ m (xi + |)] cos [k 2 2m {x 2 + b)] 

m=0 



; x g n< 0) , (38) 
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with: 

kirn = ^ , k 2m = \J ' {^f - (k\ m f , Kfc| m > , <Zk$ m > , UJ > , 

(39) 

it being understood that b 2 := {B^ ; n G N} is an unknown vector. 

6 Determination of the various unknown coeffi- 
cients by application of the boundary and con- 
tinuity conditions on Tq and I\ 

6.1 Application of the boundary and continuity conditions 
concerning the traction on Fq 



From (fT31) and (|T%|) we obtain 

rd/2 
-d/2 

,.2 



fi 1 / d X2 u 1 (x\, 0, uj) exp(— \kiixi)dxi — 

J -d/2 

/to/2 
d X2 u 2 (xi, 0, uj) exp(— \kiixi)dx\ = ; V! £ Z . (40) 
-w/2 

Introducing the appropriate field representations therein and making use of 
the orthogonality relation 

-d/2 



/ exp[i(fei n - ku)xi\dxi = d5 n i ; Vn,Z e Z , (41) 

J -d/2 



gives rise to 



oo 2 7,2 

4-^ = ^e ifc ^ 2 ^^^^/- nS in(fcl m 6) ; Vn e Z , (42) 

m=0 z " 

wherein 

^mn = / exp(±i/ci n w7/) cos(k 2 m wr])dri = 
Jo 

-e ±ifel »f {sine ± fc ln )|] + (-l) m sinc [(-k 2 lm ± fc lfl )|] } . (43) 



T 



6.2 Application of the continuity condition concerning the 
displacement on Tq 

From f|15D we obtain 

v}(x\, 0, uj) cos [kf n (xi + w/2)] dx\ — 

2 (xi,0,cj)cos [A; 2 n (xi + w/2)] dx x = ; Z = 0,1,2,.... (44) 



2 

/ - 



u 

W 
' 2 
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Introducing the appropriate field representations therein, and making use of 
the orthogonality relation 

cos [kf m (x 1 + w/2)] cos [kf n (xi + w/2)] dxi = — 5 mn 

; Vm , n = 0,1,2,.... , (45) 

gives rise to 



B ™ = TW~h\ i A n + B n] I mn e- lklnW/2 ; Vm = 0,l,2, (46) 

6.3 Application of the continuity conditions concerning the 
traction on 

From ()2U|) we obtain 

fd/2 
-d/2 



/jP / d X2 u°(xi, h, ui) exp(— ikuxi)dxi — 

J-d/2 

fd/2 

/i 1 / d X2 u 1 (xi, h,uj) exp(— \kiixi)dx\ = ; V/ G Z . (47) 



Introducing the appropriate field representations therein, and making use of 
the orthogonality relation ([41 j) . gives rise to 

- ^Aje-^ + + /i 1 *^^-^* - /i 1 ^^* = ; 

Vn G Z . (48) 

6.4 Application of the continuity condition concerning the 
displacement on I\ 

From (|19|) we obtain 

d/2 

it°(xi, /i, w) exp(— ifeija;i)da;i— 

-d/2 

d/2 

u (xi, /i,cj) exp(— ikux\)dxi = ; VIgZ. (49) 

-d/2 

Introducing the appropriate field representations therein, and making use of 
the orthogonality relation ([41 1). gives rise to 

e -ifc 2 °„h + _ ^g-i^nA _ B y k ^ h = ; Vn G Z . (50) 
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6.5 Determination of the various unknowns 

6.5.1 Elimination of to obtain a linear system of equations 
for 5°. 



After a series of substitutions, the following matrix equation is obtained 
o. 



Cl Q B Q n - D nm B° m = F n ; VnGZ, (51) 



771= — OO 

wherein: 



Cl l := cos (fcU) " i^TT sin (k\ n h) ; Vn € Z , j, i = 0, 1 , (52) 
Dutu := £ tan Ce 1 ^"^ , (53) 



and 

F n = S(u)e- [k °^ h C™8 n0 + 

^ iSj^e-^we^kl sin (fc|b) J~ J+ j( fcln _ fcl w 2 

Z=0 ^ ft 2n d 

Eq. (|51|) is a matrix equation, which, in principle, enables the determination 
of the vector b° := {B^ ; n = 0, ±1, ±2, ...}. 

6.5.2 Elimination of B® to obtain a linear system of equations for 



The procedure is again to make a series of substitutions which now leads 

a 
i 



to the linear system for B? ; V/ G N: 



B i = Pi+Y.QrniB 2 m ; V/eN, (55) 

m=0 

wherein 
and 

^ = g ^. 2e; exp (— ifc° 2 fc) exp (-ifcmf ) J f + 



71= — OO 



cos 

2f ; exn (—\ki) 

S{L0 



2ei exp {-\k\ti) exp (-ifcjf ) 



(57) 



cos [klb) C 10 

Eq. (j55fl is a matrix equation, which, in principle, enables the determination 
of the vector b 2 := {B^ ; m = 0, 1, 2, ...}. 
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7 Modal Analysis 



7.1 General considerations 

At this point, it is important to recall that the ultimate goal of this 
investigation is to predict the response of an urban site to a seismic wave. 
This response takes the form of the displacement field at various locations 
on the ground as a function of time. Thus, the field quantities of interest 
are the space-time expressions of u J ; j = 0, 1, 2. 

On account of what was written above (see (j3J) and (|24|1). the space-time 
framework diffracted field in f^o can be written as 

/OO /'OO 
dki / du> B°(k\,u) exp {i \k\x\ + k®(x2 — h) — ut\ } , 
-oo J — oo 

. (58) 

with similar types of expressions for the fields in fii and f^, as well as for the 
(incident) excitation field. These expressions possess at least four important 
features. 

The first feature, underlined in |18j and ^Hj> is that the time frame- 
work fields are expressed as integrals over the horizontal wavenumber k\ 
and angular frequency to of functions that represent, for each k\ and uj (and 
assuming there is no material attenuation in the media), either a propagat- 
ing or evanescent plane wave, the amplitude of which is a function such as 

B°(k h L0). 

A second important feature, underlined in JH]) IE] and in the com- 
panion paper, is that the amplitude functions, such as B°(ki,uj), exhibit 
resonant behavior (i.e., can become large in the presence of material losses, 
or even larger (and sometimes, infinite) in the absence of material losses), 
in the neighborhood of certain values, k* of k\ and of u, which are 
characteristic of the modes of the structure giving rise to these fields. 

The third feature, brought out in |18) . |19j . is that resonances can be 
provoked by the solicitation when: a) the frequency u of one of the spectral 
components of the latter is equal to one of the natural frequencies a;* and 
b) the horizontal wavenumber of one of the component plane waves of the 
excitation field is equal to k* . 

The fourth feature is that a mode (such as the well-known Love mode) 
corresponds to ||&;*|| larger than ||fc°||, which means that the plane wave 
associated with an excited mode is necessarily evanescent in £Iq. Conse- 
quently, to bring ||fei|| (which is a sort of momentum) up to the required 
level, requires a momentum boost, which is provided either by the incident 
field (i.e., the latter should contain evanescent wave components with hori- 
zontal wavenumbers \\ki\\ > \k°\\) and/or by the scattering structure (site) 
itself. 

In |18j . |19j . the site had horizontal, flat boundaries and interfaces and 
all the media were homogeneous, so that it could not provide the required 
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momentum boost. Moreover, it was shown in JB]> ED that if the incident 
field takes the form of a propagating (bulk) plane wave, the (Love) modes 
of the site cannot be excited, this being possible only if the incident field 
contains the required evanescent wave component, as is the situation in 
which the wave is radiated by a line source. 

Herein, the incident field takes the form of a propagating plane wave and 
the momentum boost is provided by the periodic uneveness of the surface (in 
quanta of ^), as manifested by the presence of evanescent waves in the field 
representations, so that we can expect the configuration to exhibit resonant 
behavior corresponding to the excitation of some sort of modes. 

The remainder of this section is devoted to the characterization of these 
modes and to the methods for finding the {k* ,uj*) with which they are 
associated. 



7.2 The emergence of the quasi-Love modes of the configu- 
ration from the iterative solution of the matrix equation 
for £° 

Adding and substracting the same term on the left side of the matrix 
equation (fSTj) gives 

oo 

(C™-D nn )B° = D nm B° m {l-5 mn ) + F n ; VneZ, (59) 

m=— oo 

from which we obtain 

D Em=-oo D nrn Bm (1 — 5 mn ) + F n 

B n = T^n Fi ; Vn G Z . (60) 



An iterative approach for solving this matrix system consists in computing 
successively: 

BT ] = rl0 !" ; VnGZ, (61) 

flOft) = BT ] + ^ m= -°° r Z n : VreGZ ' (62) 

and so forth. 

The Z-th order iterative approximation of the solution is thus of the form 
rO(0 ^ ._ /qoN 

wherein 

A^ 0) = F n , (64) 
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oo 

A/#>°) = F n + D nm B^ (1 - 5 mn ) , (65) 

m=— oo 

from which it becomes apparent that the solution B^ l \ to any order I of 
approximation, is expressed as a fraction, the denominator of which (not 
depending on the order of approximation), can become small for certain 
values of k\ n and u so as to make (u), and (possibly) the field in the 
substratum, large at these values. 

When this happens, a natural mode of the configuration, comprising the 
blocks, the soft layer and the hard half substratum, is excited, this tak- 
ing the form of a resonance with respect to B® , i.e., with respect to the 
field in the substratum. As B® is related to A\ and B\ via (|151) - (|5n|) . the 
structural resonance also manifests itself in the layer for the same k\ n and u. 

Remark 

The matrix equation (|51j) can be written as 

(C-D)b° = f (66) 

wherein: the matrix C 10 has components C^ m = C*°<5 nm , D has components 
D nrn , and the vectors b° and f have components B® and F n respectively. The 
modes of the configuration are obtained by turning off the excitation |38| . 
embodied in the vector f . Thus, the non-trivial solution of the homogeneous 
matrix equation ()66|) is the solution of (the dispersion relation) 

det (C 10 - D) = (67) 

wherein det(M) signifies the determinant of the matrix M. 

The equations 

■D n = C™-D nn = ; n = 0,±l,.... , (68) 

associated with a singularity in the iterative procedure (|63|) . also correspond 
to an approximation of the dispersion equation of the modes (|67[) when the 
off-diagonal elements of the matrix C 10 — D are small compared to the 
diagonal elements. Such a situation does not necessarily prevail, but it 
is nevertheless useful to obtain a first idea of the natural frequencies of 
the modes from the simple relations (|68|) rather than from the much more 
complicated relation (|6*7)) . 

As shown in the companion paper (section 6.2), and in |18[ I19j . it is 
impossible to excite a Love mode in a configuration without blocks consisting 
of a soft layer overlying a hard halfspace when the incident wave is a plane 
bulk wave. This case corresponds to 6^0. 
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Let us return to the denominator of the expression of B„, which takes 



the form: 

V n = - D nn = C£"> - £ ^T«t tan (k^b) . (69) 

l=Q V k 2n d 

Remark 

Cf = , (70) 
is the dispersion relation for ordinary Love modes. 

Remark 

V n = 0, (71) 

is then the dispersion relation of what we term quasi-Love modes which are 
generally different from ordinary Love modes. 

Remark 

When b — > 0, the dispersion relation for quasi-Love modes becomes the dis- 
persion relation T> n = = for ordinary Love modes. 

Remark 

fA^ k^ w 

For small P 1 , 2 1 m , , the quasi-Love modes are a small perturbation of ordi- 

V k 2n d 

nary Love modes. 

7.3 The emergence of the quasi displacement-free base block 
modes and quasi-Cutler modes of the configuration from 
iterative solutions of the linear system of equations for 

R 2(0 

7.3.1 Approximate dispersion relations arising from the first type 
of iterative scheme 

Let us consider (|55|). which can be re- written as: 

00 

Bf (1 - Qu) = + Qml (1 - 8ml) B m \ V! £ N . (72) 
m=0 

A (first type of) iterative procedure for solving this linear set of equations 
leads to: 

1 - Qu 

i?, 2(1) M = Bf\u) + E "=° Qml (1 ~ 6ml) B ™ 0) ; 1 = 1,2 (74) 

l - Qu 
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This procedure signifies that becomes large when 1 — Qu is small, and 

that this occurs at all orders p of approximation. The fact that B^' be- 
comes inordinately large is associated with the excitation of a natural mode 
of the configuration. The equations 1 — Qu = ; I G N are the approximate 
dispersion relations of the l-th natural modes (I G N) of the configuration. 

Remark 

We say that these dispersion relations are approximate in nature because 
they are obtained by neglecting the off-diagonal terms in the matrix equa- 
tion ()72|) . This may not be legitimate, but it is nevertheless useful to get an 
idea of the true natural frequencies by examination of the solutions of the 
approximate dispersion relations 1 — Q mm = ; m G N. 
Let us therefore examine the latter in more detail: 

00 1,2 2 f-YOl 

^ m (fci, w) := cot [kl m b) - £ i €m ^^^I- n I+ = ; m G N , (75) 

n=— oo n2^ n 

which shows that the modes of the configuration result from the interaction 
of the fields in two substructures: the superstructure (i.e., the blocks above 
the ground), associated with the term 

Tx m (ki,u)) = T\ m = cot {kl m b) , (76) 

and the substructure (i.e., the soft layer plus the hard half space below the 
ground) associated with the term 

n=-oo a K n2^ °n 

Each of these two substructures possesses its own modes, i.e., arising from 
Ji m = 0, for the superstructure, and Tim = 0, for the substructure, but the 
modes of the complete structure are neither the modes of the superstructure 
nor those of the substructure, since they are defined by 

F lm (h,uj)- Tton{ki,u) =0 ; meN , (78) 

which again emphasizes the fact that the modes of the complete structure 
result from the interaction of the modes of the component structures. 

In order to obtain the natural frequencies of the complete structure, we 
first analyze the natural frequencies of each substructure, and assume that 
all the media are non-dissipative (i.e. elastic). 

The solutions of the (approximate) dispersion relations for the super- 
structure are: 



/, 9 ,\ , 9 , (2n + l)vr 
cot (kl m b) = k 2 2m b = i —!— 




(2n + l)ir\ 2 fm-K\ , . 

[ —^) +(— ) ; n,m = 0,l,2,..., (79) 
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which are the natural frequencies of vibration of a block with displacement- 
free base (i.e., at these natural frequencies, v?\ x _ Q vanishes on the base 
segment of the block). 

Next consider the dispersion relations for the geophysical (sub)structure 
J~2m{k\,to) = ; m £ N. As pointed out in the companion paper, the sum 
in this relation can be split into three parts corresponding to: i) propaga- 
tive waves in both the substratum and the layer, ii) evanescent waves in the 
substratum and propagative waves the layer, iii) evanescent waves in both 
the substratum and layer. Only the second part can lead to a vanishing 
denominator, and also to the satisfaction of the dispersion relation of Love 
modes. 

Remark 

2 

For small jp, the quasi displacement-free base block modes are a small per- 
turbation of the displacement-free base block modes. 

Remark 

For small ^f, the quasi displacement-free base block modes are a small per- 
turbation of the displacement-free base block modes. 

Remark 

We notice that the approximate dispersion relations for the configuration 
involving an infinite set of equispaced identical blocks are similar to the dis- 
persion relations obtained (in the companion paper) for a small number of 
blocks, in that they betray the existence of a combination of quasi-Love and 
quasi displacement-free base block modes, which are small perturbations 

of Love and displacement-free base block modes respectively for small 

M 

and/or small ^. 

7.3.2 Approximate dispersion relations resulting from a second 
type of iterative scheme 

Let Q denote the matrix of components Q m u I the identity matrix, and 
b and p the vectors of components Bf and Pi respectively. 

The system of linear equations (|55[) can be written as the matrix equation 
(for the determination of the unknown vector b) 

(I - Q) b = p . (80) 

The modes of the configuration are obtained by turning off the excitation 
|37j . embodied in p. The non-trivial solutions of l|8fljl are then obtained from 

det (I - Q) = . (81) 
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An iterative (partition) procedure for solving this equation, which is differ- 
ent from the preceding one, and is particularly appropriate if the off-diagonal 
elements of the matrix are small (but non neglected) compared to the diago- 
nal elements, is first to consider the matrix to have one row and one column, 
i.e., 

1 - Qoo = , (82) 
then to consider it to have two rows and two columns, 

det(jr Qo ° ~ Qm ) =(l-Qoo)(l-Qn)-QoiQio = 0, (83) 

and so forth. 

We shall not go into the details of these various approximate expressions 
of the dispersion relation because they become extremely involved beyond 
(|82|). However, it is of historical and didactic interest to examine (|82|) in 
detail. This is done in the next section. 



7.3.3 Solution of the zeroth-order dispersion relation arising in 
the two types of iterative schemes: the Cutler mode 

We rewrite the lowest-order approximation of the dispersion relation (|80[) 
(equivalent to for n = 0) in the form first given in |36j : 

oo -7,2 2 /"»01 

Vo = l~ £ tan (k 2 b) I^Xn^to = • (84) 

n=— oo 2nr^ u 

However: 

T~ = sine 2 I 

2 



hJL = sinc2 ( k in^r J , (85) 



so that 

„ w k 2 u? , „ N ^-^ k° 9 / w\ C® 1 , , 

^0 = 1- i-T^TT tan (fc 2 6) £ ^sinc 2 (fc ln -) ^ = . (86) 



d fcV 



We now consider the cases (first studied in [TT1 I3T1 131 123 1 H\ ?H\ ) in 
which the layer is filled with the same material M° as that of the substratum 
(actually, in most of the cited publications, M 2 was also taken equal to M°, 
but this is not done here, for the moment at least). Thus, /U 1 = /jP and 
k 1 = k°. In addition, we recall that M° is non-dissipative, so that fjP and 
A; are real. We shall also suppose, to simplify matters, that fj? and k 2 are 
real (i.e., M 2 is lossless). 
Then 

C = C™ = e W (ikLh) = exp^/i) =>■ ^ = 1 , (87) 
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W kZfJ, 2 , ,,\ \^ k° 2/7 w \ 

2>0 = 1 - 1 7P^o tan ^ 6 )z^yr smc (^7)=° 



so that 

If, in addition, we assume that M 2 = M°, then 

V Q = 1 - tan (fc°6) ^ ^sinc" = , (89) 

with 7] = 2. This dispersion relation is identical to the one obtained in [Q 
(wherein 77 is not given a definite value). 

Let us stick to (|88(l for the moment and make the substitutions 

6:= 1 ' a:= F' K:= W = V ' ^ = d' (90) 



therein, so that 

A; 

^2n 



£> = 1 - l^Ktan (fc 2 6) V To-smc 2 (fe ln — j = . (91) 



Suppose that there exists an integer N for which one of the terms in the 
series is very large compared to the rest of the series. We thus write 



Vq = 1 — i(5Ktan (k 2 b) 

with 



^2N 



W 



sine ( kiN— I + R\ 



, (92) 



Rq 



S IcT sinc2 ) • ( 93 ) 

ne ( Z _{AT}) 2n 



Let fci := fcijv and &2 := fc^iV- It follows that 

7T WW WW 

k ln = k! + 2(n-N)- fc Xn - = A*- + (n- JV)?r- = A*- + (n- JV)tt<5 . 

(94) 



Now, for the n = N term to be large, at the very least 

« 1 . (95) 



u w 
2 



and 

k%N ~ ■ (96) 
The first of these two conditions is synonymous with 

sinc^l) «1. (97) 
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A corollary is that 



sin (ki n — j = sin ( fa — - ) cos ((n — N)tt5) + 



cos (fa~) sin ((n - iV)mf) « sin ((n - N)nS) . (98) 
We shall also assume that 

8^1 , (99) 



so that 



whence 



sin (kir^j ~ sin ((n - iV» ss => R ss , (100) 

1 - Was tan (£: 2 &) ^-sinc 2 (fa-) = , (101) 

v k,2 v 2 / 

which is the dispersion relation obtained by Rotman |31| . 
It is more precise to take account of (|9"T|l so that 

k° 

1 - i5 k tan (k 2 b) — = , (102) 

which, when k = 1 (i.e., the case M 2 = M°), constitutes the dispersion 
relation first obtained by Cutler jllj . We shall now examine this relation in 
detail. 

Since all the parameters in the Cutler dispersion relation are real, the 
latter has no solution unless &2 = y (k ) 2 — (ki) 2 is imaginary, i.e., k2 = 
ill 1 / (k 1 ) 2 — (/co) 2 ||) which occurs if \ki\ > k°. If we recall that imaginary ki 
corresponds to an evanescent (i.e., surface) wave, then we can say that the 
Cutler mode is associated with the excitation of a surface wave. 

We now inquire as to the conditions under which the Cutler mode can 
be excited. The first step is to find the values of (k\,uj) which are solutions 
of 

k° 

l-5Ktan(k 2 b) -— = ; \ki\ > k° . (103) 
|| fall 

Another point of view is to consider fa to be the wavenumber (of a surface 
wave) to which is associated the phase velocity ci such that 

fa = — . (104) 
ci 

Then it is easy to obtain (from (|103|0 



ci = c° ./ ? . (105) 

l/ l+[^tan(f)] 2 



Remark 
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This relation, first published in [2^, shows that even if M 2 is non-dispersive 
(recall that it was assumed, from the start, that M° is non-dispersive), then 
the phase velocity of the Cutler mode is dispersive, i.e., c\ = c\(oS) (which, 
of course, is the reason why one speaks of a dispersion relation in connection 
with a (e.g., Cutler) mode). 

Remark 

The Cutler mode corresponds to a slow (surface) wave with repect to the 
bulk plane waves in Qq, since c\ < c°. 

Remark 

c\{uj) is a periodic function of u, since c\{uj + ^§-) = c\(uj) ; VI 6 N. 
Remark 

Ci {uj) = c° for a£ = Itt ; V/ G N. 
Remark 

ci(fa;) = for ^ = (21 + l)f ; VieN. 
.Remark 

The phase velocity of the Cutler mode is all the closer to the phase velocity 
of bulk waves in M° for all to, the smaller is k. On the contrary, for a Cutler 
mode with phase velocity very different from that of bulk waves in wc 
must have a large contrast between the material properties of M° and M 2 . 

Let us now inquire as to the means of actually exciting a Cutler mode 
with an incident plane bulk wave. At first, this seems impossible (for the 
same reason it is not possible to excite a Love mode with an incident plane 
bulk wave). But we must not forget that the field in f^o is composed not 
only of diffracted plane bulk waves, but also of diffracted evanescent waves, 
the possibility of these waves to exist being due to the uneven (at present, 
periodically uneven) geometry of the stress-free surface at the site. 

The discussion concerning the Cutler mode began with the assumptions: 
i) the term in the expression of the diffracted field in $1° corresponding to 
the iV-th order diffracted plane wave dominates all the other terms, and ii) 
this diffracted wave is an evanescent wave, i.e., k^M 1S imaginary. In the 
dispersion equation context, k± is a variable that has no particular connec- 
tion with the solicitation. When the site is solicited by a plane bulk wave, 
then k\N = k\ + %-iV, with k\ = k° smO 1 the factor directly related to the 
solicitation (9 % the angle of incidence). Thus, for the iV-th order evanescent 
diffracted wave to be excited, we must have 

2ir 

kx = k 1N ^kx = k\ + —N , (106) 
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with N such that fcj^y = i|| y^fciAr) 2 — (/c°) 2 ||. This so-called coupling re- 
lation, i.e.. ()l(J6j) . translates the fact that the periodic topography adds the 
momentum ^-N necessary to convert the incident bulk wave into an evanes- 
cent wave (whose phase velocity is smaller than that of the bulk wave, and 
whose horizontal wavenumber is therefore larger than the wavenumber k° 
of the incident bulk wave) . 

In order for a resonance to occur in the iV-th order mode, the frequency 
of one of the components of the spectrum of the excitation must be equal to 
a natural frequency of the mode. Although this is a necessary condition, it is 
not a sufficient condition, because we must also have ||£;i|| = ||£;i7v|| > \\k°\\. 

A remark is in order concerning what happens when Rq is not neglected 
in the expression of T>q. The subset, in this remainder term, involving 
the n for which ft2 n is imaginary (i.e., corresponding to evanescent waves), 
will modify somewhat the (real) solutions of what formerly constituted the 
Cutler dispersion relation, whereas the subset involving the n for which k2 n 
is real (i.e., corresponding to propagative waves) adds an imaginary part 
to the (real) solutions of what formerly constituted the Cutler dispersion 
relation. Thus, the evanescent Cutler wave becomes a leaky wave, i.e., a 
wave with complex k\. 

7.3.4 Solution of the zeroth-order dispersion relation when M 1 7^ 
M°: the quasi-Cutler mode 

The dispersion relation (|84|) has been studied in |36| and is a subset of 
the dispersion relations analyzed in sect. 17.21 Not much more, other than 
what is revealed by a numerical analysis, can be added to the text in sects. 
l7~2l and 17X31 due to the complexity of (fHlj). 

To make a long story short, one finds that the presence of the layer 
transforms the Cutler mode into a quasi Cutler mode which is all the closer 
to a Cutler mode the smaller is the layer thickness h and/or the closer the 
material parameters of the layer are to those of the substratum. If, on the 
other hand, h is not very small, and/or the material parameters of the layer 
are very different from those of the substratum, then the quasi-Cutler mode 
becomes an entity entirely different from that of the Cutler mode; in fact, 
it ressembles a Love mode, so that it is better to represent the phenomena 
in terms of quasi-Love modes (as in sect. 17. 2j) than in terms of quasi-Cutler 
modes. 

8 Computation of the fields w°, u l and u 2 

8.1 Computation of u° 

The quasi-modal coefficients -B^(w), Vm 6 N are obtained by employing 
the partition procedure (i.e., reducing the infinite-order matrix in (|80|) to 
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a M x M matrix and the vectors to M-tuple vectors, solving the finite- 
order matrix equation so obtained, and increasing M until convergence is 
obtained of the successive M-th order approximate solutions). Once the 
B^lo), Vm G N are computed in this manner, the field in the block domain 
O2 is obtained via (j38j) . This field vanishes on the ground at the frequencies 
of occurrence of the displacement-free base modes of the block. 

8.2 Computation of u 1 

Let us next consider the field in the layer. Combining ([42)1 . (|46j) . (|48|) 
and (|5U|) leads, via (|3"5|) . to: 

1 2S(co) exp (i (k\xi — kn 2 h)) cos (^02^2) 
u (x,w) = -p + 



00 00 



n=-oom=0 r" n2 n 

,01.0 



cos (fc n2 (x 2 - h)) + i 1 " 2 sin (k n2 (x 2 - h)) 



■ (107) 



8.3 Computation of u 2 



Let us finally consider the field in the substratum. Combining (|42l) , ()46l) , 
(|i8|) and (JSOl) leads, via (EU, to: 

u°(x, cj) = u*(x, w) + 

cos (k\ 2 h) + i^mr sin (kl 2 ti) 
S(u) exp (i (fcixa + k° 02 (x 2 - 2/i))) + 

00 

^2 iexp (i (fcnixi + k° n2 (x 2 - /1))) X 



ra=— 00 



E^^f^^l^exp^l) . (108) 



8.4 Comments on the fields u l and u 2 

Expressions (|107j) and (|108j) indicate that both displacement fields 
u 1 (x, uj) and u 2 (x, w) are composed of: i) the field obtained in the absence 
of the blocks and induced in the layer or substratum by the incident plane 
wave, ii) the field induced by the presence of the blocks, which appears as 
a field radiated by an infinite number of identical source distributions (each 
one related to a given block). In particular, each of these induced sources 
takes the form of a ribbon source of width w located at the base segment of a 



30 



block when it is related to the zeroth-order quasi-mode (see the companion 
paper) . 

When a mode is excited (i.e., at a resonance frequency), one or several of 
the can become large, in which case it is possible for the fields to become 
large at resonance. This will be demonstrated in the numerical examples 
which follow. 

9 Numerical results for the seismic response in 
two idealized cities 

9.1 Preliminaries 

The numerical results are obtained in two manners: i) by the Mode- 
Matching (MM) method (described in the previous section) as it applies to 
configurations consisting of an infinite number of equally-spaced, equally- 
sized rectangular blocks, and ii) by the Finite-Element (FE) method (briefly 
described in the companion paper, and in more detail in JfHEl) as it applies 
to configurations with a large (but finite) number of equally-spaced, equally- 
sized rectangular blocks. 

The discussions concerning the numerical aspects of the dispersion rela- 
tions will be based on material stemming from the MM. For the purpose of 
the analysis, and to allow for an easier comparison with the results exposed 
in the companion paper, we re- write (the lowest-order approximation of the 
dispersion relation) (|H4|) in the form 



9.2 Ten and an infinite number of blocks in a Nice-like site 
9.2.1 Parameters 

The sources of earthquakes in the city of Nice (France) are usually deep 
and located almost vertically below the Nice site (a few kilometers laterally 
from the city). Modeling the solicitation of the site by a normally-incident 
(i.e., 9 l = 0, which leads to k\ = 0) plane incident wave is thus realistic |32| . 
The central frequency of the Ricker pulse associated with the solicitation is 
chosen to be vq = 2Hz. 

With Nb the number of blocks, we consider here both the cases Nf, = 10 
and Nb = oo. 

The parameters of the underground of our idealized Nice urban site are: 
p° = 2200 Kg/m 3 , c°=1000 m/s, Q° = oo, p 1 = 1800 Kg/m 3 , cJ e/ =200 m/s, 
Q 1 = 25, with the soft layer thickness being h = 50 m. 
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v (Hz) v (Hz) 

Figure 6: Indications concerning the solution of the dispersion relation T\ — 
J~2 = for the Nice-like city. In the left panel: the solid and dashed curves 
depict 3i(.Fi) and ^(^2) respectively versus frequency [y in Hz) . In the 
right panel: the solid and dashed curves depict ^(^i) and ^(^"2) respectively 
versus frequency. 



The Haskell eigenfrequencies, which are usually close to the quasi-Love 
mode frequencies, are then v^ ASK = 2 ™ +1 1^ =lHz, 3Hz, 5Hz,.... 

The blocks are chosen to be b =30m high, w =10m wide, and their 
center-to-center spacing is d =50m. Their material constants are chosen to 
be: p 2 = 250 Kg/m 3 , c 2 e/ =240 m/s, Q 2 = 10. 

Thus, the displacement-free based block eigenfrequencies (solutions of 
T 2 = 0), v® F = (2m + l)c 2 46, are v® F = 2Hz, GHz... and the quasi-Cutler 
mode natural frequencies (which are specific to the periodic nature of the 
site, with characteristic dimension d), are obtained from J 7 = T\ — J-i = 0. 

9.2.2 Infinite number of blocks 

Fig. El gives an indication of the natural frequencies of the modes of 
the global configuration. Recall that such a natural frequency is a solution 
of T\ — Ti = 0, requiring that ^(^i) = ^(.T^), at the least. This occurs 
at v « 2.5Hz, 5.1Hz,.... in the frequency range of the figure. The atten- 
uation associated with a particular mode (at a frequency v*) is related to 

9(^(1/*)) -9f(^ ? (^)). 

One observes in the figure that the quasi-Love modes, which are expected 
to occur near 1Hz, 3Hz, 5Hz,... are either not excited or are strongly at- 
tenuated, while what appears to be the quasi-displacement free base block 
mode at i>Q DFB ps 2.5 is excited with a relatively-low attenuation. On 
the contrary, what appears to be the quasi stress-free base block mode at 
v qdfb ra g \ff z ^ \ s associated with a large attenuation and should therefore 
have little effect on the global response of the site. 

The notable features of this response are that it is dominated by the 
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Figure 7: 2tt times the spectrum (left panel) and time history (right panel) 
of the total displacement at the center of the top segment of a block, which 
is the leftmost one of a set of 10 blocks (solid curves, obtained by the FE 
method), or any one of an infinite set of blocks (dashed curves, obtained by 
the MM method ) in a Nice-like site. 



quasi-displacement-free block modes, and that the influence of the periodic 
nature of the distribution of blocks, which manifests itself by the quasi- 
Cutler modes, probably appears at frequencies higher than those in the 
range of the figure. 

9.2.3 Comparison of responses for Nb = 10 and Nb = oo 

We now compare both the spectra and time histories for an infinite 
number of identical blocks with those of a finite number (i.e., 10) of the 
same blocks. The point of observation is at the center of the top segment 
of a block. The block under examination is the leftmost one in fig. Q and 
the fifth one in fig. |HJ Since the angle of incidence is 0, both the spectra 
and time histories of response do not change from one block to another in 
the Nb = oo structure. We note the close similarity of both the frequency 
and time domain responses relative to the Nb = 10 and Nb = oo structures, 
even for the leftmost block of the finite configuration. This means that for a 
Nice-like site, the model involving an infinite number of identical blocks can 
used to model and/or to analyze the response phenomena in such a city. 

Other than this, as concerns the response spectra, the left-most peak at 
~ 1.1 Hz and the rightmost peak at ~ 3.2 Hz seem to be related to the 
excitation of the Haskell pseudo-modes (which is another way of saying that 
the corresponding quasi-Love modes are weakly excited, since these peaks 
would exist even if the buildings had nearly zero height). The large peak 
in the response spectra at ~ 2.5Hz is associated with the excitation of the 
quasi-displacement free base block mode, as expected from examination of 
the solutions of the dispersion relation (see fig. EJ). The lengthening of 
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Figure 8: 2tt times the spectrum (left panel) and time history (right panel) 
of the total displacement at the center of the top segment of a block, which is 
the fifth one of a set of 10 blocks (solid curves, obtained by the FE method), 
or any one of an infinite set of blocks (dashed curves, obtained by the MM 
method ) in a Nice-like site. 



the duration (« 8 sec) of the response, with respect to the duration of the 
incident pulse of ~ 1 sec, is mainly due to the excitation of this mode, but 
since the quality factor of the resonance of this mode is relatively small, 
the lengthening of the duration is relatively- modest. Thus, it seems that 
both the amplification of the amplitude, and the lengthening of duration, of 
seismic motion in a periodic or quasi-periodic portion of Nice, cannot take 
catastrophic proportions. Insofar as this conclusion concerns the amplitude 
of ground motion, it is in agreement with the conclusions in |32| . 

9.3 Finite and infinite number of blocks in a Mexico City-like 
site 

9.3.1 Parameters 

The sources of the major earthquakes in Mexico City have usually been 
shallow and located in the subduction zone off the Pacific coast, approxi- 
mately 350km from Mexico City), so that modeling the solicitation of this 
city by a plane incident bulk wave is not realistic Nevertheless, we 

assume such a (normally-incident) plane bulk wave solicitation, notably to 
enable an easy quantitative comparison between the finite and infinite city 
responses and qualitative comparison with previous studies [38 l ITU l 132} 14^ IK] . 
The central frequency of the Packer pulse associated with the solicitation is 
chosen to be uq = 0.5Hz. 

The underground of the city is characterized by: p° = 2000 kg/m 3 , 
c°=600 m/s, Q° = oo, p 1 = 1300 kg/m 3 , c^ e/ =60 m/s, Q 1 = 30, with the 
soft layer thickness being h = 50 m. The Haskell frequencies are 0.3Hz, 
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0.9Hz, 1.5Hz,.... 

The blocks are 50m high, 30m in width, and their center-to-center spac- 
ing are successively chosen to be 65m, 150m and 300m. The material con- 
stants of the blocks are: p 2 = 325 Kg/m 3 , c;; e/ =100m/s, Q 2 = 100. The 
natural frequencies of the displacement-free base block modes are v^f B = 
0.5Hz, 1.5Hz,..., and the quasi-Cutler mode natural frequencies (which are 
specific to the periodic nature of the site, with characteristic dimension d), 
are obtained from J- = 0. 

9.3.2 Dispersion characteristics of the modes for an infinite num- 
ber of blocks 

Fig. Ogives indications concerning the solutions of the dispersion rela- 
tion of the global configuration (i.e., blocks plus underground) for center-to- 
center spacings d = 65m, d = 150m and d = 300m. For a center-to-center 
spacing d = 65m, the influence of the periodic nature of the city (embod- 
ied in the parameter d) appears essentially at a high frequency outside the 
spectral bandwidth of the solicitation, while for d = 150m and d = 300m, 
this influence can make itself felt in the spectral range of the solicitation 
through the quasi-Cutler modes. Quasi-Love and quasi-displacement free 
base block modes should be excited at 0.3Hz and « 0.5Hz respectively, 
both associated with a low attenuation. The stress-free base block mode is 
excited at ~ 1Hz, but is highly-attenuated for all three d. 

Quasi-Cutler modes are clearly excited with a low attenuation. In par- 
ticular, for a relatively large center-to-center spacing (i.e. d = 300), the 
spectral density of this type of mode is large. The fundamental quasi-Love 
natural frequency is probably very close to the fundamental quasi-Cutler 
natural frequency. 

9.3.3 Seismic response for an infinite number of blocks with pe- 
riod d = 65m compared to that of no blocks 

In fig. El] we compare the spectra and time histories in the presence and 
absence of the d = Q5 set of blocks. 

The first resonance peak in the spectra is shifted to a lower frequency 
and is of higher amplitude when the blocks are present. This is the in- 
dication of the excitation of the fundamental quasi-Love mode responsible 
for what was termed the soil- structure interaction in the companion paper. 
The peak relative to the excitation of the perturbed displacement-free base 
block mode has a relatively high amplitude, due to both the spectrum of 
the incident wave and to the presence of the other blocks. Effectively, the 
resonance frequency does not correspond to that of the displacement-free 
base block mode (at which frequency the response is nil at the center of the 
base segment of the block), and the particular shape of the spectrum at the 
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Figure 9: Indications concerning the solution of the dispersion relation T\ — 
J~2 = for the Mexico City-like urban site. In the left panels: the solid and 
dashed curves depict ^(J^i) and ^(^-"2) respectively versus frequency (u in 
Hz) . In the right panels: the solid and dashed curves depict ^(Fi) and 
^s{J--2) respectively versus frequency. The top, middle and bottom panels 
are related to d = 65m, d = 150m and d = 300m respectively. 
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Figure 10: Comparison of 2ir times the spectrum (left panels) and time 
history (right panels) of the total displacement on the ground in absence of 
blocks (solid curves) with the displacement at three locations in the presence 
of blocks (dashed curves): i) on the ground, at a location halfway between 
two adjacent blocks (bottom panels), ii) at the center of the base segment of 
a block (middle panels) , iii) at the center of the summit segment of a block 
(top panels). 
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center of the base segment of the block is characteristic of the excitation of a 
multi-displacement-free base block mode and not to the excitation of a quasi 
displacement-free base block mode (as shown in the companion paper). This 
coupled mode also takes into account the so-called structure-soil-structure 
interaction as (see the companion paper) was suggested in sect. 17.31 

A splitting of the first peak appears in the ground (between the blocks) 
response. The second of the split peaks has a low quality factor, which fact 
means the quasi-absence of beatings in the time history of response. 

The temporal displacements are of higher amplitude and duration in the 
presence of blocks than in the absence of blocks, particularly at the center of 
the top segment of the block. These results are evocative of those obtained 
in the companion paper for two blocks. 
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9.3.4 Comparison of the seismic responses for one, two and an 
infinite number of blocks with center-to-center separations 

d = 65m 

In fig. 1111 we compare the spectra and time histories of response for 
one-block, two-block, and infinite-block configurations at the center of the 
top segment of one of the blocks. This figure shows that the effect of in- 
creasing the number of blocks is essentially to increase the height of the 
second (i.e., higher-frequency) resonance peak at the expense of the first 
resonance peak, and thus to introduce stronger high-frequency oscillations 
in the temporal response. As the quality factor of the second resonance peak 
increases with Nb, the durations also increase with iV&. The combined effect 
of the increased quality factors and higher frequency oscillations is increased 
cumulative motion of the block. 

In fig. 1121 we compare the spectra and time histories of seismic response 
for one-block, two-block, and infinite-block configurations at the center of 
the bottom segment of one of the blocks. Since the second resonance now is 
synonymous with a minimum of response, the aforementioned effects are not 
produced at the base of the block. In fact, with increasing Nb, we actually 
observe a decrease of the height of the first resonance peak, whose effect is 
to decrease the amplitude and duration of motion in the time histories. 
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Figure 11: Comparison of 27r times the spectrum (left panels) and time his- 
tory (right panels) of the total displacement on the ground in absence of 
blocks (solid curves) with the displacement at the center of the top segment 
of a 50m x 30m block (dashed curves) for: i) one such blocks (top panels), ii) 
two such blocks with center-to-center spacing d = 65m (middle panels) , and 
iii) an infinite number of such blocks with center-to-center spacing d = 65m 
(bottom panels). The one- and two-block results are taken from the com- 
panion paper and apply to the response to the wave emitted by a deep line 
source situated at (0m, 3000m) whose amplitude, at ground level is different 
from that of the plane wave soliciting the infinite block structure. 
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Figure 12: Comparison of 2tt times the spectrum (left panels) and time 
history (right panels) of the total displacement on the ground in absence 
of blocks (solid curves) with the displacement at the center of the bottom 
segment of a 50m x 30m block (dashed curves) for: i) one such blocks (top 
panels), ii) two such blocks with center-to-center spacing d = 65m (middle 
panels), and hi) an infinite number of such blocks with center-to-center 
spacing d = 65m (bottom panels). The one- and two-block results are taken 
from the companion paper and apply to the response to the wave emitted by 
a deep line source situated at (0m, 3000m) whose amplitude, at ground level 
is different from that of the plane wave soliciting the infinite block structure. 
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9.3.5 Comparison of the seismic responses for ten, twenty, fourty 
block configurations with that of a configuration having an 
infinite number of blocks for center-to-center separations 

d = 65m 

In fig. El we compare the spectra and time histories of seismic response 
at the center of the top segment of a centrally-located block in configurations 
with 10, 20, 40 and an infinite number of 50m x 30m blocks separated by 
d = 65m. On the whole, these displacement responses are all the same, in 
both the frequency and time domains, marked by relatively-long duration 
duration (~ 2 min), and large maximum and cumulative amplitudes. How- 
ever, for the finite values of iV&, there appears some splitting of the Nb = oo 
low frequency resonance peak which gives rise to beatings in addition to 
those due to the presence of the two main resonance peaks. 

In fig. El w e compare the spectra and time histories of seismic response 
at the center of the bottom segment of a central block in configurations 
with 10, 20, 40 and an infinite number of 50m x 30m blocks separated by 
d = 65m. The splittings referred-to in the previous lines are now more 
apparent and give rise to more pronounced beatings, especially for the 10 
and 20 block configurations. These beatings are absent for the Nb = oo city. 
The durations are quite long for the 10 and 20 block configurations, and, on 
the whole, the signals are evocative of what has been often observed during 
earthquakes in certain districts of Mexico City. 

In fig. El we compare the spectra and time histories of seismic response 
at the midpoint on the ground between two adjacent centrally-located blocks 
in configurations with 10, 20, 40 and an infinite number of 50m x 30m blocks 
separated by d = 65m. Again, we observe two main resonance peaks for all 
Nb, giving rise to characteristic beatings, to which are added other beatings 
due to splittings of the first resonance peak, especially noticeable for finite 
Nb- This response is again quite evocative of ground response observed in 
the midst of certain districts of Mexico City during many earthquakes that 
have affected this city jHj. 

The presence of the additional low-frequency peaks in this set of figures 
is linked either to the (finite) number of blocks considered and/or to the 
total width W of the finite configuration (for N = 10, 20, 40, W is equal 
to 650m, 1300, 2600m respectively). These peaks cannot be accounted-for 
in the dispersion relations written above, since they result from an analysis 
of configurations with an infinite number of blocks (and for which W = oo). 
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Figure 13: Comparison of 2-rr times the spectra (left panels) and time histo- 
ries (right panels) of the total displacement at the center of the top segment 
of a 50m x 30m block of a configuration with an infinite number of blocks 
(solid curves) with the displacement at the same location of a central block 
(dashed curves) of configurations having: i) ten such blocks (top panels), ii) 
twenty such blocks (middle panels), and iii) fourty such blocks. In all cases, 
the center-to-center spacing is d = 65m. 
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Figure 14: Comparison of 2tt times the spectra (left panels) and time his- 
tories (right panels) of the total displacement at the center of the bottom 
segment of a 50m x 30m block of a configuration with an infinite number of 
blocks (solid curves) with the displacement at the same location of a cen- 
tral block (dashed curves) of configurations having: i) ten such blocks (top 
panels), ii) twenty such blocks (middle panels), and iii) fourty such blocks. 
In all cases, the center-to-center spacing is d = 65m. 
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Figure 15: Comparison of 2ir times the spectra (left panels) and time histo- 
ries (right panels) of the total displacement at the midpoint on the ground 
between adjacent 50m x 30m blocks of a configuration with an infinite number 
of blocks (solid curves) with the displacement at the same location (dashed 
curves) of a central region of configurations having: i) ten such blocks (top 
panels), ii) twenty such blocks (middle panels), and iii) fourty such blocks. 
In all cases, the center-to-center spacing is d = 65m. 
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9.3.6 Illustration of the spatial variability of response in a con- 
figuration of ten blocks for center-to-center separations d = 

65m 

Fig. 1161 depicts snapshots of the displacement field in the entire configu- 
ration containing ten blocks whose center-to-center distance is d = 65m. It 
can be noticed that: i) the motions are still quite strong in portions of the 
configuration some 50s after the arrival of the initial pulse and ii) that the 
motion is quite variable spatially- and temporally-speaking, iii) the spatial 
variability extends even to within a given block. 

The spatial variability of response is one of the characteristics of the 
motion that has often been recorded within Mexico City |14[ 18]. 
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Figure 16: Snapshots at instants t = 25s (top panel), t = 32.5s (middle 
panel), t = 50s (bottom panel), of the total displacement field for 10 identical 
50m x 30m blocks whose center-to-center spacing is d = 65m. 
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9.3.7 Comparison of responses of the configuration without blocks 
to the one with iVj, = oo blocks for d = 300m 

In fig. I17( we compare the spectra and time histories in the presence and 
in the absence of the blocks. When the blocks are present, their number is 
infinite and their center-to-center distance is 300m. 

The excitation of the multi displacement-free base block mode cannot be 
clearly distinguished because of the excitation of quasi-Cutler modes whose 
existence is related to the periodic nature of the block distribution. Note 
should be taken of the fact that now (i.e., for d = 300m) these modes are 
visible (with a large quality factor) within the bandwidth of the solicitation, 
whereas they were invisible for the d = 65m periodic structure. As previ- 
ously, the soil-structure interaction appears as a resonance peak associated 
with the excitation of the fundamental quasi-Love mode. 

These features show up at all three locations of the configuration with 
blocks. They manifest themselves in the time domain by: i) a larger duration 
(multiplied by ~ 4 on the top segment of the blocks with respect to its value 
in the absence of the blocks), ii) a larger peak amplitude (at the top of the 
blocks) and larger cumulative motion, and iii) pronounced beatings at all 
locations due to the periodic nature of the configuration. 

These features are once again evocative of those which have been ob- 
served during earthquakes in certain districts Mexico City. However, the 
excitation of the quasi-Cutler mode, which is strongly-linked to the quasi- 
periodic or periodic nature of a district of the city (these districts actually 
exist, as seen in fig. can, at best, explain only part of the features of 
response in this city, since the latter is not usually solicited by a (normally- 
incident) plane wave. 

It has been shown in the companion paper, and in |181ll9j that the correct 
solicitation of a configuration for the study of the Michoacan earthquake 
(and many other earthquakes affecting Mexico City) is a cylindrical wave 
radiated by a shallow, laterally-distant source that gives rise to Love waves 
after traveling within the crust from the hypocenter to the city. In the 
same studies, it was shown that this type of solicitation is another cause for 
the large duration, large amplitude, and beatings of Mexico City response. 
Unfortunately, source wave solicitation cannot be treated by the quasi-modal 
analysis of periodic structures. 



48 



50 r 



e 40 



\ 30 



20 



27T||U Z (X,V)|| 
27T||U 1 (X,V)|| 




50 r 



E 40 



1 30 



20 



27t||u' ! (x,v)|| 
- 2 It ||u 1 (x,v)|| 




50 r 

40 

30 



infinite number of 50mx30m blocks 
without blocks 



JL 20 



10 




0.2 0.4 0.6 0.8 ,1 ,1.2 1.4 1.6 1.8 2 
v(Hz) 



0.2 0.4 0.6 0.8 ,1, .1.2 1.4 1.6 1.8 2 

v(Hz) 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

v(Hz) 



6 

E 4 

i- 2 

IT. 

x -2 
^ -4 
-6 



6 
4 

"E 2 
E, 

s; 
- 3 -2 
-4 
-6 



u 2 (x,t) 
u 1 (x,t) 




20 40 60 80 100 120 140 160 
t(s) 



u 2 (x,t) 
u 1 (x,t) 




20 40 60 80 100 120 140 160 
t(s) 



infinite number of 50mx30m blocks 
without blocks 




20 40 60 m 100 120 140 160 

t(s) 



Figure 17: Comparison of 2tt times the spectrum (left panels) and time his- 
tory (right panels) of the total displacement on the ground in the absence 
of blocks (solid curves) with the displacements (dashed curves) at the cen- 
ter of the top segment (top panels), center of the bottom segment (middle 
panels) and midpoint on the ground between adjacent 50m x 30m blocks, 
of a Nb = oo configuration in which d = 300m. 
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9.3.8 Comparison of the MM and FEM results for a finite set of 
blocks whose separation is d =300m 

We now compare, in figs. El El and EDI the results computed by the 
mode-matching technique (accounting for the fundamental quasi-mode of 
the blocks), for an infinite number of equally-spaced identical blocks, with 
those we obtain by our finite-element method, for a finite number = 10 
of equally-spaced identical blocks. In this case, B = {1,2, 10}, with the 
block numbers 1 and 10 being located at the left and right edges respectively 
of the finite configuration. 

The width of the finite set of blocks W is 3000m and therefore closer 
(than in the previous case) to the infinite width implicit in the MM theory. 
This may be the reason why the FE results agree quite well with the MM 
results within the city. The two computational modes match less well at a 
point on the ground outside the city (fig. l2*U|) . as one would expect. 

One observes that the response is spatially variable, of long duration 
(attaining at some points ~ 2 min), with high maximum and cumulative 
amplitudes, and characterized by beating features, which is quite evocative 
of the features of many sismograms of Mexico City earthquakes (HI I14| I33[ 
I15j . This suggests that rather largely-spaced blocks or buildings are more 
apt than closely-spaced blocks or buildings to induce the large-scale features 
(particularly the very long durations and beatings) observed in this city. 
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Figure 18: Comparison of 2n times the spectrum (left panels) and time 
history (right panels) of the total displacement as calculated by the MM 
method for an infinite number of 50m x 30m blocks separated by d = 300m 
(solid curves) with those obtained by the FEM method in the central portion 
of a configuration of 10 identical 50m x 30m blocks separated by d = 300m 
(dotted curves) : i) at the center of the top segment of block number 5 (top 
panels), ii) at the center of the base segment of block number 5, and iii) on 
the ground between block numbers 5 and 6. 
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Figure 19: Comparison of 2n times the spectrum (left panels) and time 
history (right panels) of the total displacement as calculated by the MM 
method for an infinite number of 50m x 30m blocks separated by d = 300m 
(solid curves) with those obtained by the FEM method in the left-hand edge 
portion of a configuration of 10 identical 50m x 30m blocks separated by 
d = 300m (dotted curves): i) at the center of the top segment of block 
number 1 (top panels) , ii) at the center of the base segment of block number 
1, and iii) on the ground between block numbers 1 and 2. 
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Figure 20: Comparison of 2n times the spectrum (left panels) and time 
history (right panels) of the total displacement as calculated by the MM 
method for an infinite number of 50m x 30m blocks separated by d = 300m 
(solid curves) with those obtained by the FEM method in the left-hand edge 
portion of a configuration of 10 identical 50m x 30m blocks separated by 
d = 300m (dotted curves) at a point on the ground 150m to the left of the 
block 1. 
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Figure 21: Snapshots at various instants, t = 25s, t = 32.5s, t = 50s, of the 
total displacement field for 10 identical equally-spaced 50m x 30m relatively 
distant blocks configuration, solicited by a normal incident plane wave. The 
center-to-center spacing is d = 300m. 



9.3.9 Illustration of the spatial variability of response in a con- 
figuration of ten blocks separated by d = 300m 

Fig. contains three snapshots, at instants t=25s, 32.5s and 50s, of 
the displacement field in the entire configuration of 10 blocks with center- 
to-center separation of 300m. Once again, one notes both the large- and 
small-scale spatio-temporal variability of response in the city site. 

10 Conclusions 

The response, to a plane wave initially propagating in the substratum, 
of a finite set of identical, equi-spaced blocks, each block modeling one, or 
a group of buildings, in welded contact with a soft layer overlying a hard 
half space, was investigated in a theoretical manner via the mode matching 
technique. 

The capacity of this technique to account for the complex phenomena 
provoked by the presence of blocks on the ground was demonstrated by 
comparison of the numerical results to which it leads to those obtained by 
a finite element method. 

It was shown that the presence of the blocks induces a modification of 



54 



the phenomena that are produced by the configuration without blocks, or 
of a configuration of closed blocks disconnected from the geophysical half- 
space. In particular, the blocks modify the dispersion relation of what, in the 
absence of the blocks, constitutes the Love modes. Moreover, the periodic 
nature of the urban site is responsible for the existence of another vibrational 
mode which is a variant of the Cutler mode encountered in electromagnetic 
waveguide structures. 

The dispersion relation for the periodic configurations with an infinite 
number of blocks is very complex, but an approximation of this relation 
lends itself to a fairly-explicit theoretical analysis, inspired notably by the 
method first adopted in the electromagnetic wave community. The general 
features of this dispersion relation, revealed by the theoretical analysis (and 
manifested by the existence of several types of vibrational modes), were 
shown to actually exist by means of the numerical study. 

The excitation of the vibrational modes was then studied in the par- 
ticular case of city districts with 10 and an infinite number of identical, 
periodically-disposed blocks. A common feature of the influence of the 
blocks is the excitation of the quasi-Love mode, which occurs even for so- 
licitation by a plane wave (recall that, for this type of solicitation, it is not 
possible to excite ordinary Love modes in a flat ground (i.e., no blocks) /soft 
horizontal layer /hard substratum configuration JB])- The trace of quasi- 
Love mode excitation in the frequency domain was shown to be a shift 
to lower frequency and an increase of the amplitude of the first (lowest- 
frequency) peak of the response. 

The change of the phenomena provoked by plane wave solicitation, from 
a configuration without blocks (for which there exist only bulk waves in 
the geophysical structure), to one with blocks (for which there exist quasi- 
Love modes characterized by a field in the substratum that is predominantly 
a surface wave in the substratum) is a manifestation of the so-called soil- 
structure interaction. 

Multi displacement-free base block modes were shown to be excited in 
all the configurations and to correspond to a coupled mode. 

The modifications of the response due to the presence of blocks separated 
by 65m were found to be rather modest in the Nice-like site. 

The modifications of the frequency-domain response (and less so of the 
time-domain response) were found to be fairly substantial for structures in- 
volving blocks separated by 65m in the Mexico City-like site. In particular, 
10 or 20 blocks separated by 65m, were shown to give rise to anomalous 
features (amplifications of peak and cumulative motion, large durations and 
beatings) that are even closer to those observed in Mexico City than config- 
urations with a larger number (40, oo) of blocks. 

All the anomalous features found for the Mexico City-like site with 10 or 
an infinite number of blocks were found to be closer to the actually-observed 
anomalous features observed during earthquakes in Mexico City when the 



55 



separation between blocks is 300m rather than 65m. This was found to 
be due to the fact that the structure with the larger period enables the 
quasi-Cutler mode to make itself felt within the range of frequencies of the 
source. 

The theoretical findings and numerical results of the present study can- 
not account for all the anomalous phenomena observed in many of the earth- 
quakes in Mexico City (notably the exceptionally-long durations) for the 
obvious reasons that the model adopted herein, i.e., SH plane wave solicita- 
tion, 2D, periodic geometries, simple underground (horizontal homogeneous 
soft layer of infinite lateral extent overlying a homogeneous lossless hard 
substratum), simple homogenized building blocks, linear soil behavior, etc., 
is incomplete, and in some respects, rather far removed from reality. Never- 
theless, the results of our study indicate that it is possible that the excitation 
of vibrational modes, whose structure is closely-related to those described 
herein, was responsible for at least part of the large-scale, anomalous me- 
chanical effects that have caused so much damage in past earthquakes in 
urban areas. 

The most important finding of this work, which substantiates those ob- 
tained in [HH1IS11E1, is that the presence of groups of buildings (i.e., city 
blocks) can modify substantially the seismic motion in a city. Moreover, 
provided the blocks are arranged quasi-periodically with a sufficiently large 
period d, the seismic motion is of longer duration, and of higher cumulative 
(sometimes peak) amplitude on the ground (and, of course, in the buildings) 
than when the built structures are not present. 

Therefore, it is advisable to integrate (as is starting to be done [T3] ) 
the presence, composition and layout, of buildings, together with, and to 
the same extent as, the features of the underground structure and compo- 
sition, into the large-scale computer codes that are being employed to 
predict the level and durations of shaking in highly-populated, economically- 
important, earthquake-prone areas. 
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